Aging in Models of Non-linear Diffusion 



o 

CD 

Q 

0^ 



> 

00 

O 



X3 
O 

o 



X 
S3 



Daniel A. Stariolo 
Dipartimento di Fisica, Universita di Roma I "La Sapienza" 

and 

Istituto Nazionale di Fisica Nucleare, Sezione di Roma I 
Piazzale Aldo Mow 2, 00185 Rome, Raly 
stariolo@chimera . romal . inf n . it 

We show that for a class of problems described by the non- linear diffusion equation ^< 
an exact calculation of the two time autocorrelation function gives C{t,t') = f{t — t')g{t'), {t > t') 
exhibiting normal and anomalous diffusions, as well as aging effects, depending on the values of ^ 
and V. We also discuss the form in which the Fluctuation- Dissipation Theorem is violated in this 
type of systems. Finally, we argue that in this kind of models, aging may be consequence of the non 
conservation of the "total mass" . 



D- 



In a wide variety of physical systems where some kind 
of diffusion takes place it can be observed that the mean- 
squared displacement scales with time as < (t) > cx 
with a depending on the physical problem in question. 
a = 1 corresponds to the so called normal diffusion (the 
simple random walk), of which a complete statistical de- 
scription can be obtained, for example, from the solu- 
tion of the well known diffusion equation ■^(j){x,t) = 

D-^4>{x, t), where t) is the probability that the dif- 
fusing particle be at position x at time t provided it was 
at the origin x — a,t t — and D is the diffusion con- 
stant. If a ^ 1 the diffusion is called anomalous with 
a < 1 corresponding to suhdiffusion and a > 1 to su- 
perdiffusion 0| . Anomalous diffusion can be consequence, 
for example, of some kind of disorder in the system [|l],D, 
or more generically, of longe range correlations in space- 
time. The computation of the propagator (f>{x,t), which 
contains all the spatio-temporal information of the sys- 
tem, is in general a difficult task. Without knowning the 
exact propagator for all times, its long time form can be 
calculated in some cases using techniques like renormal- 
ization group and scaling arguments |^ . 

Recently Tsallis et.al. Q have obtained the exact so- 
lution of the non-linear Fokker-Planck equation 
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{x,tr = -^{F{xmx,t)r}+Df 



dx^ 



(1) 



where (/x, G TZ^, D > is a diffusion constant, 
F{x) — —dV{x)/dx is an external force associated with 
the potential V{x) and (x,t) is 1-1-1 space-time. They 
have found the solution for a drift of the form F[x) = 
k\ — k2X with ki and k2 constants. This equation re- 
covers the standard diffusion or Fokker-Planck equation 
when ^ = V = 1. Other values of (/i, v) represent interest- 
ing physical systems as well: in the case with F{x) = 
(a purely diffusive problem), for fi = 1 and arbitrary 



V, Eq.(|]) is known as the Porous Medium Equation and 
models many non-equilibrium systems in fluid dynamics 
particle diffusion in magnetic fields and gas dy- 
namics , depending on the value of v. The information- 
theoretic aspects of this case have been studied by Plas- 
tino and Plastino |^ and recently the interplay between 
dynamic and thermodynamic aspects have been studied 
in [|lO| for the general case of diffusion in N dimensions. 
The case ^jl — 1 and v = ?> has been studied by Spohn ||] 
and describes a solid-on-solid model of surface growth. 

Restricting to the situation without drift (for the gen- 
eral case see Ref. Q), the solution for the propagator 
4>q{x,t) can be written as 



<t>qix,t) 



[l-/3(f)(l-g)[x-XM (t)] 
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(2) 



with q = I + ji — v and XM{t) = X]\i(0) is the mean 
position, which for a situation without drift is constant 
and equal to the initial position. This solution is closed 
by the relations satisfied by /3(t) and Zq(t), namely: 
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and 



[^,(0)]^ 



z,{t) 



2ly{ly + ^l)D(3{0)[Zg{0)]^^' 



(3) 



(4) 



A static form of Eq.(H) with /3(t) = 1/T (inverse temper- 
ature) and Zq{t) — Zq{T) (partition function) has been 
obtained from a maximum entropy principle in the con- 
text of a generalized thermstatistics |l^ , and succesfully 
applied for explaining, among many other problems, the 
thermodynamic foundations of Levy anomalous diffusion 
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|^,|3[. We will see in the following that the above so- 
lution presents a very rich dynamical behaviour charac- 
terized in general by anomalous diffusion and, for certain 
values of /i and by aging phenomena, the long term 
memory effects observed and nowadays extensively stud- 
ied in amorphous polymers [Q and spin glasses ||l5| . Let 
us consider the two times autocorrelation function 

C{t,t')^ {y{t)x(t')) (5) 

dxdyxy(l)q{x,0,t')(l)q{y,x,t- t') , (6) 

in which t' < t and where (j}q{u, v, z — z') is the probability 
that the particle was at position u at time z provided it 
was at position v at time z' . From Eq. (H) we obtain 



CM^Kq Zqit-t')Zqit'mt~t')nm?^^ 



(7) 



with Kq a constant that only depends on q. Now consid- 
ering the regime in which t — t' ^ oo and also t,t' — > oo, 
from Eqs.(|), (§ and (|) 



C{t, t')^A[B{t- t')] ^ [B t'] ^ 



(8) 



where 



A 



and 



/32(0) Z4a'(0) 
(9) 

(10) 



This result presents a variety of interesting charac- 
teristics. First we note that for ^ — v = 1 we ob- 
tain the well known result for normal diffusion C(<, t') = 
2 Dmin{t,t'). For the case of the porous medium equa- 
tion, i.e. fj, = I, the asymptotic correlation simplifies 
to 



C{t,t') =A[Bt']~ 



(11) 



In this case the long time behaviour depends only on the 
minimum time (as in normal diffusion) but the diffusion 
is anomalous with exponent 2/{l + ly). When v > 1 the 
behaviour is subdiffusive and for < 1 it is superdiffusive. 
This qualitative change can be conveniently observed in 
the shape of the propagator at a fixed time t* as shown 
in Figure 1. In the subdiffusive regime, characteristic of 
the porous medium equation, the propagator presents a 
"wave front" which expands with time as 
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v//3(0)(i.- 1)^,(0) 



(12) 



with Zq{t) given by Eq.(^. As i ^ cxo it flattens as 
\xwf\ oc The "wavefront" changes to the exponen- 

tial decay of the normal gaussian diffusion when v = 1 
and then to superdiffusion characterized by power law 
tails in the propagator which decays for x ±cx) and t 
fixed as 



(x) oc \x\ 



(13) 



Although the previous discussion considered /it = 1, it is 
valid in the general case in which the different regimes are 
more conveniently characterized by the parameter q = 
— i/(g<l,l and > 1 corresponds to subdiffusive, 
normal and superdiffusive behaviours respectively), and 
the relation defines the possible choices for and v. 

For [I ^ \ the autocorrelation function depends ex- 
plicitly on both times, on all time scales, a feature char- 
acteristic of systems with long term memory. These ef- 
fects of aging are common in disordered systems (like, 
e.g. spin glasses) where the time correlations present 
particular scaling forms (see |l5|] and refernces therein). 
Aging effects have also been studied in models without 
explicit disorder as scalar fields, the XY model and spin- 
odal decompositio n and also in random walks in dis- 
ordered media | |T7| , |l8| . In order to compare the present 
scalings with those found in spin glasses it is convenient 
to normalize the correlation function in order to have 
C{t,t) = 1. This normalization is natural in magnetic 
systems but it is not so for an arbitrary random walk. In 
our problem a suitable definition may be 



Cn{t, t') 



C{t,t') 



y/C{t,t)C{t',t') 

With this definition we obtain 



, M-1 /'t'\ 2(f'+.-) 



(14) 



(15) 



when t^t' 00 and t — t' — > oo. This normalized cor- 
relation presents the commonly seen scaling with t' /t. 
The breaking of time translation invariance can be clearly 
seen in the so called "displacement" function 

B{t,t')^{W)-x{t')f) (16) 
^C{t,t) + C{t\t')-2C{t,t'). (17) 

By definition B{t, t) = 0. Writing t = t-t': 
B{t' + T,t') (xt'TT^ -^(^t' + t)tt^ _2T7S^i'w (18) 

We see that in general B{t' + t, t') depends on both t' and 
T and consequently time translation invariance is broken. 
When — 1 and t ^ t' 



B{t'+T,t') « T~ 



(19) 



The displacement depends only on the time difference t 
and time translation invariance is recovered. 
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These results also reflect in the behaviour of a suitable 
generalization of the Fluctuation-Dissipation Theorem 
(FDT) for systems that may never attain equilibrium. 
This generalized form, firstly introduced for studying the 
off-equilibrium dynamics of spin glasses states that 

R{t,t')^l30{t^t')X{t,t')^^^^. (20) 

where R{t, t') is the response at time t to an external 
force h applied at time t' 



R{t,t') ^ 



S{x{t)) 



5h{t') 



(21) 



h=0 



P is the inverse temperature and 6{z) is the step function. 
The function X{t,t') measures the departure from FDT; 
if FDT is satisfied then X{t, t') = 1. In general FDT will 
be violated by a system that never reaches equilibrium as 
in the models we are considering here, but it is neverthe- 
less instructive to analyze how the function X behaves. 
If we apply an external perturbation h at time t' , the 
form of the propagator remains the same as in Eq.(^, 
but now XM{i) satisfies the differential equation 



dt 



XM{t)^h6{t-t'), 



(22) 



whose solution is xmit) = xm{0) + h9{t — t'). According 
to Eq.(^l|) the response function is R{t,t') = d{t — t'). 
Now, for the case of the usual random walk (/x = = 1), 
an explicit calculation from Eq.(^ gives for the correla- 
tion function C{t,t') — 2Dt'. Noting that in the stan- 
dard Fokker-Planck Equation f3 — > 1/D, we obtain from 
Eq.(|0|) that X{t,t') = 1/2. Consequently, for the sim- 
ple random walk, FDT is violated with a factor X that 
is a constant 0. In the case of the Porous Medium 
Equation {fJ. — 1, arbitrary v), a similar analysis shows 
that X{t,t') — f{t'), i.e., the correction factor is a func- 
tion of t' only. In the third case, for arbitrary fj, and v 
the function X depends explicitly on both times t and 
t' . From this analysis it is clear that, although all three 
cases violate FDT (as expected), the dynamics becomes 
more complex as we go from the simple random walk, 
to the Porous Medium Equation and finally to arbitrary 
ifj-jiy), and this is reflected in the structure of the function 
X{t,t'). 

Finally and from other point of view, we argue that 
the absence (presence) of aging phenomena is a direct 
consequence of the "conservation of mass" (or not) in 
the system. For our model it can be verified that M 



dx 4>q(x, t) 



dx4>q{x,Q), (23) 



and consequently the norm of the distribution or "total 
mass" is independent of time only if /i = 1 . As suggested 
by Goldendfeld for the case of the normal diffusion, the 



conservation of mass is directly connected with the fact 
that asymptotically the system loses memory of its his- 
tory and becomes unable to distinguish between different 
initial conditions. The situation is completely different if 
the mass is not conserved, as for example, in the modi- 
fied porous medium equation and in the problems we 
consider here when /i 7^ 1. Although we are not able to 
give a rigorous proof of this assertion, the different situa- 
tions studied here are in agreement with it. Consequently 
it would be a direct connection between conservation of 
the mass and aging in the class of systems considered: 
for our model equation in the systems with /i = 1 the 
mass is conserved and the autocorrelation function does 
not present aging behaviour. But if 7^ 1 there is no 
conservation of the mass and aging effects are seen. 

Concluding, a number of interesting problems that 
can be modeled by non-linear diffusion equations can 
be solved exactly and the diffusion presents very dif- 
ferent characteristics depending on the degree of non- 
linearity. An exact calculation of the two time correla- 
tion functions shows that, besides the expected anoma- 
lous diffusion, some systems may exhibit aging effects as 
found in many other disordered systems. The aging sat- 
isfies particular scaling forms, depending on the problem 
considered, which in principle could be compared with 
experimental results. The different dynamical scenar- 
ios can be characterized also studying the (violation of) 
Fluctuation-Dissipation Theorem. These aging effects, in 
systems modeled by partial differential equations, seem 
to be strongly related to the conservation of the "total 
mass". It would be interesting to test the validity (or 
not) of this hypothesis in other models and also to study 
the extensions of these results for dimensions higher than 
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FIG. 1. Qualitative shape of the propagator (j)q{x,t) for 
fixed t. From top to bottom q = —1, 0, 1, 2, 3 



[1] Bouchaud J P and Georges A 1990 Phys. Rep. 195, 127. 

[2] Zumofen G, Klafter J and Blumen A "Models for Anoma- 
lous Diffusion" , in Disorder Effects on Relaxational Pro- 
cesses, Richert R and Blumen A (Eds.), (Springer- Verlag, 
Berlin, 1994). 

[3] Goldenfeld N Lectures on Phase Transitions and the 
Renormalization Group, (Addison- Wesley Pub. Co., se- 
ries Frontiers in Physics, 1992). 

[4] Tsallis C and Bukman D J 1996 Phys. Rev. E 54, R2197. 

[5] Aronson D G "The Porous Medium Equation", in Non- 
Linear Diffusion Problems, Fasano A and Primicerio 
M (Eds.), Lecture Notes in Mathematics vol. 1224 
(Springer- Verlag, Berlin, 1986) and references therein. 

[6] Berryman J G 1977 J. Math. Phys.18, 2108. 

[7] Zeldovich Y B and Raiser Y P Physics of Shock Waves 
and High-temperature Hydrodynamic Phenomena, vol. II, 
(Academic Press, New York, 1966). 

[8] Plastino A R and Plastino A 1995 Physica A111, 347. 

[9] Spohn H 1993 J. Phys. I France S, 69. 
[10] Compte A and Jou D 1996 J. Phys. A2Q, 4321. 
[11] Tsallis C 1988 J.Stat.Phys.52, 479 

Curado E M F and Tsallis C 1991 J. Phys. A24, L69 
Corrigenda 1991 J. Phys. ^24,3187; 1992 J. Phys. 
A25,1019. 

[12] Tsallis C, de Souza A M C and Maynard R "Derivation of 
Levy- Type Anomalous Superdiffusion from Generalized 
Statistical Mechanics", in Levy Flights and Related Top- 



4 



